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Abstract. We establish analogues of the Hawking and Penrose singularity theorems based on (a) aver- 
■ aged energy conditions with exponential damping; (b) conditions on local stress-energy averages inspired 

by the Quantum Energy Inequalities satisfied by a number of quantum field theories. As particular ap- 
plications, we establish singularity theorems for the Einstein equations coupled to a classical scalar field, 
which violates the strong energy condition, and the nonminimally coupled scalar field, which also violates 
the null energy condition. 

bJQ, 

1 Introduction 

The singularity theorems proved by Penrose and Hawking in the mid-1960's [TJ [H |3] represented 
a major advance over previous work, which was restricted to situations of high symmetry and 
particular matter models. At the technical level, this advance was achieved by significant devel- 
opment of the geometrical apparatus; at the physical level, however, the key idea was to replace 
particular matter models by generic energy conditions, which drive geodesic congruences to focal 
points. While many classical matter models obey the classical energy conditions, such as perfect 
fluids and the electromagnetic field, there are exceptions. For example, the Klein-Gordon field 
with nonzero mass fails to satisfy the strong energy condition (SEC; the hypothesis for the Hawk- 
• • , ing and Hawking-Penrose theorems if the cosmological constant vanishes) while the nonminimally 

coupled Klein-Gordon field also violates the null energy condition (NEC; the hypothesis of the 
Penrose theorem). See [I] for a critical assessment of the status of classical energy conditions 
and [5] for a review of singularity theorems and commentary on energy conditions in that context. 

One of the main limitations of the energy conditions is that they are incompatible with quantum 
field theory, as has long been known [6] . Examples of states with locally negative expectation values 
for the energy density are easily constructed [7] and in a large class of quantum field theories the 
energy density at any given point is unbounded from below as the quantum state varies (see [8] 
and references therein). Moreover, recent results in two-dimensional conformal field theory show 
that individual measurements of weighted spacetime averages of energy density in the vacuum 
state yield negative values with probabilities approaching 90% in some cases [5]. 

For these reasons a number of authors have investigated singularity theorems under weakened 
energy conditions. Examples include [THl EJ HH E3] , in which various averaged energy conditions 
are considered, particularly the averaged null energy condition (ANEC), i.e., 

WVdA>0, 

where 7 is a complete (or sometimes future- or past-complete) null geodesic in an affine param- 
eterisation. [Via the Einstein equations, these conditions translate immediately to conditions on 
the Ricci tensor.] The definition of the integral requires some care; it is usually interpreted as a 
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liminf of the integrals over finite regions (as in [T^]) or regularised by inserting a mollifying factor 
/(A/Ao) 2 in the integrand (for / in a suitable function class) and taking the liminf as Ao — > oo 
[13] . On the other hand, a different approach (see 7.21-22 of [Tl] and [TS]) is to suppose a lower 
bound on R a bt a t b for all unit timelike vectors and sufficient initial contraction for focussing to 
occur. 

In this paper, we will establish singularity theorems under conditions that combine ideas from 
these two approaches but also extend them to allow for situations in which the integrand becomes 
unboundedly negative. Roughly speaking, we replace the ANEC integral by the requirement that 
for a future-complete null geodesic 7 : [0, 00) — > M, 



should be finite for some c > (with the integral understood as above) and dominate the initial 
expansion of suitable future-complete null congruences in spacetime M. This leads to a general- 
ization of the Penrose theorem, which is stated precisely as Theorem 15.21 a similar generalization 
of Hawking's theorem is given as Theorem 15.11 In Section [51 we apply these results to prove a 
Hawking singularity theorem for the massive minimally coupled Klein-Gordon field (even though 
it can violate the SEC), at least provided that the field amplitude is bounded or, more generally, 
of exponential type. Under similar conditions on the amplitude of the field and its derivatives, we 
also prove a Penrose singularity theorem for the nonminimally coupled Klein-Gordon field (even 
though it can violate the NEC). Our results are based on criteria, derived in Section [21 for the 
nonexistence of global solutions to Riccati equations that generalize arguments of |16j by inclusion 
of the exponential as an integrating factor (one could equally consider other functions in this con- 
text). This idea is used to generalize arguments of |13) in Section [3] and the relevant singularity 
theorems are proved in Section [5] 

In Section [4[ we consider how exponentially damped energy conditions of the type required 
can be derived from bounds on local energy averages similar to those satisfied by quantum fields, 
known as Quantum Energy Inequalities (QEIs) or Quantum Inequalities (QIs). Inequalities of 
this type, which express limitations on the magnitude and duration of violations of the classical 
energy conditions, were first mooted by Ford |17) and have since been established in a number of 
quantum field theory models (see [18] for reviews). In four-dimensional Minkowski space, for 
example, the energy density of a free real scalar field obeys the bound 



for any smooth compactly supported real- valued function / and any Hadamard state ip [HI .20 O 
Here 7 is any complete timelike geodesic in a proper time parameterisation. Under rescaling of /, 
the bound scales as the inverse fourth power of the sampling timescale r; this is consistent with 
the unboundedness below of the energy density at points (in the limit t — > 0) and also gives an 
averaged weak energy condition in the limit r — > 00. Similar results are known for other free fields. 

Although analogous results hold in curved spacetimes, the bounds are more complicated in 
form. For example, the bound given in |22j involves local Fourier transforms of distributions 
formed from the first few terms in the expansion of the Hadamard parametrix, and is valid for 
averaging within suitable domains that typically have compact closure. However, it is expected 
on physical grounds [23 (and borne out in concrete examples) that averages over timescales 
smaller than local curvature length scales are constrained by bounds taking the same form as in 
Minkowski space. In the bounds of [22] . for example, the first term in the Hadamard expansion will 
dominate the others on sufficiently small scales, and becomes well- approximated by the expression 
in Minkowski space. This motivates the study of energy conditions in which local averages of the 
energy density or other stress-energy quantities are bounded from below in terms of L 2 -norms 
of derivatives of the averaging function, perhaps with constraints on the support or after the 

3 Units with h = c = G = 1 are used throughout. Our geometrical conventions are those of 1211 . 
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subtraction of some reference function (as in the so-called difference QEIs) . In Section 0] we will 
derive (variants of) the exponentially damped energy conditions from local energy conditions of 
this type. This prepares the way for a more detailed study, which would have to quantify the 
scales over which the QEI bounds can be replaced by these simpler estimates. We also emphasize 
that, while our local conditions are motivated by QEIs, they do not coincide directly with them; 
see Section [6] for discussion on this point. 



2 Criteria for nonexistence of global solutions to the Riccati 
equation 

We begin with a slight modification of an argument given in [161 Lemma 3]. 
Lemma 2.1. Consider the initial value problem 

z 2 

z = \-p, 

q 

z(0) = z (2.1) 

where q{t) and p(t) are continuous on [0,oo) 7 and q(t) > on [0,oo). // 

f°° dt . f T 

/ — —r = +oo and liminf / p(t) dt > —zn 

Jo q(t) T^+ooJ Q 

then (|2.ip has no solution on [0, oo). 

Proof. Suppose there is a solution z(t) on [0, oo). By hypothesis, there exists t\ > such that 

p(t') dt' > -Z Q 





for all t € [ti, oo). Integrating the differential equation, 

f ^r d *'+ f t p(t')dt' + z > f'^-dt' (2.2) 

Jo q(t') Jo Jo q(t ) 

for t > t\. Introducing R(t) = J Q z(t') 2 /q(t') dt' , we see that R is nonnegative and obeys the 
differential inequality 

q q 

for t>t\. Accordingly, we have R(t) > for all t>t\. Fixing any ti > t\, we now have 

1 1 1 /"* R , f l dt . . 

>1^-1^= -5S*t> ~ ( 2 - 3 ) 



R(t 2 ) ~ R(t 2 ) R(t) J t2 R 2 4 q 

for t > t 2 - However, the right-hand side is unbounded as t — > oo and we obtain a contradiction. □ 
Here, 'no solution' means that we have z(t) — > +oo as t — > t~ < oo. Note that 

z(t) >z +[ p(t') dt' 







for all t for which the solution exists, so z{t) — > — oo cannot occur at finite time. (It is easily shown 
that divergence is the only way in which the solution can break down at finite times.) 
As a digression from our main development, we note the following consequence. 
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Corollary 2.2. //, in Lemma \2.1[ the integral condition on p is replaced by 

f T 

z + inf / p(t) dt = a > , 



then (|2.1I) has no solution on [0, t], where r is the unique solution to 

<- T dt' _ 2 
o ?(*') a ' 

In particular, this applies if the negative part p~(t) = min{0, p(t)} is integrable and zq+ p^(t)dt - 
a > 0. 

Proof. Following the proof of the Lemma, we may take t\ = and deduce from the hypothesis 
and Eq. (12.21) that z(t) > a for all t £ [0, oo) and hence that 



R{t 2 ) > a 



2 , t2 dt' 



lit') 

for any t 2 > 0. The inequality (12.31) now asserts that 

1 f f t2 dt' \ / ? dt' 

^7 > 



o q(t')J \J t2 q(t' 

for all < ti < t and the intermediate value theorem allows us to find £2 such that both factors 
on the right-hand side are equal to | J * q(t')~ 1 dt' , establishing the required result. Finally, if p_ 

is integrable, it is clear that J Q T p(t)dt > J Q °° p_(t)dt for all T. □ 

Remark. In the standard singularity theorems, q(t) = n — 1 for timelike geodesic congruences and 
q(t) = n — 2 for null geodesic congruences, where n is the spacetime dimension. The upper bound 
on the parameter time until blow-up is therefore r = 2(n—l)/a or 2(n — 2)/a respectively. Tighter 
bounds can be obtained in this case, by a more involved argument 24 . 

We now weaken the conditions on p\ for simplicity we also fix q to be a constant. 
Lemma 2.3. Consider the initial value problem 

z 2 

z = hr. 

s 

2(0) = z (2.4) 
where r{t) is continuous on [0, oo), and s > is constant. If there exists c > such that 

z Q -- + lim inf / e' 2ct/s r(t) dt > 
then (|2.4p has no solution on [0, 00). 

Proof. Again, assume the contrary. Then y(t) = (z(t) — c)e~ 2ct / s solves 



y 



y 



2 



+ e- 2ct / s (r(t)+c 2 /s) 



se -2ct/s 

2/(0) = z - c (2.5) 



on [0, 00). This equation is of the form ((271]), with q(t) = se~ 2ct/s and pit) = e~ 2ct / s (r(r) + c 2 /s). 
Clearly, J °° dt/q{t) = 00; as we also have 

lim inf / e- 2ct/s (r(t) + c 2 1 s) dt > liminf / e - 2ct/s r(t) dt + liminf / e~ 2ct f s c 2 1 s dt 

T^+ooJ T^+ooJ T->+ooJ Q 

= - + lim inf / e- 2ct/s r(t)dt 
> c- z Q = -y(0) 
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by the hypothesis, Lemma I2TT1 entails that (|2.5[) has no solution on [0, oo) and we have obtained a 
contradiction. □ 
For example, suppose that r(t) > —Ae Bt for A,B > 0. Then for c > Bs/2 we have 



' liminf C e- 2ct/s r{t)dt>- C -- 



2 t^+oo J w " 2 2c-5s 

the right-hand side of which has a maximum of — (\/As + Bs/4) at c = VAs + Bs/2. Accordingly, 
(12. 4[) has no solution on [0, oo) if 

z a > VAl: + — . (2.6) 
4 

As an explicit example, consider r(t) = — e 2t , s = 1 (i.e., A = 1, B = 2). The above inequality 
establishes nonexistence of global solutions for zq > 3/2; in fact, the equation can be solved in 
terms of modified Bessel functions as 

z(t) = j t log (Me^K^l) + KaWhO)] - [/o(e*)#o(l) - ifo(e%(l)]z ) 

and it is easy to see that z(t) has a singularity at finite time if and only if 

1^)^(1) + K (e t )I 1 (l) _ KM 
20 > l>o I (e*)K (l) - K (et)I (l) K (l) ' 

which is consistent with the estimate provided by Lemma 12.31 
We also have the following consequence. 

Corollary 2.4. If, in Lemma \2.3[ the integral condition on r is replaced by 

r T 

zo-c+inf / e- 2ct/s (r(t)+c 2 /s)dt = a>0, 

T>0j 

then (|2.4)l has no solution on [0, r], where 

T = Yc log ^ + 4c / a ) • 

In particular, this applies if the negative part r_(t) = min{0,r(t)} is such that r_(t)e~ 2ct / s is 
integrable and zq — c + j °° e~ 2ct '< /l 'r ■-(t)dt = a > 0. 

Proof. Following the proof of Lemma 12.31 the hypothesis implies that 

y(0) + inf / pit) dt = a>0 

T>oJ 

and we may apply Corollary 12.21 to find that there is no solution to Eq. (|2.4j) on [0,t], where r 
solves 

2 f T dt' e 2cT l s - 1 



a J q(t>) 2c 

proving the first statement. The second statement follows from the elementary bound e~ 2ct l s {r(t)- 
c 2 /s)dt > J Q °° e- 2ct ' s r^{t)dt for all T E [0, oo). □ 
For cases in which a 3> c, the bound is approximately 2s/ a, which is the same as in the remark 
following Cor. [221 
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3 A generalization of the Wald—Yurt sever argument 

Wald and Yurtsever [T3] consider hypotheses on weighted averages of the forcing term in the 
Riccati eqution (see the Lemma in that reference). Here, we generalize this further, after the 
pattern of Lemma 12.31 above. Fix a real- valued compactly supported function g £ (R) with 
g(t) = 1 on [0, 1], and g nonincreasing on R + . 

Lemma 3.1. Suppose that a solution to (|2.4[) exists on (all of) t £ [0, oo). If r(t) is nonconstant 
then 



limsup / e- 2ct/s r{t)g{t/T) 2 dt < - - z 

T— >0O JO * 



for all c > 0. If r(t) = r is constant, the same conclusion holds provided c ^ zq or c 2 /s + r ^ 0. 
Remarks. 

1. As the liminf of any function is always less than its limsup we may rather trivially replace 
limsup by liminf in the statement of this result. 

2. The c — case essentially corresponds to the result in [13] . 

Proof. If a solution to (|2.4|) exists for all t, then we know that (|2.5|) has a solution on [0, oo). 
Rearranging, and integrating against g(t/r) 2 , 

c 2 f°° 



h(r)~ e- 2ct ' s r{t)g{t/r) 2 dt= y( t )g(t/r) 2 dt- - e^ ' s g{t I r) 2 dt - R c (r) , (3.1) 

JO JO s J0 

where, as before, y(t) = (z(t) — c)e~ 2ct / s obeys Eq. (|2.5|) and 

1 f°° 

Ro{t) = - V (t) 2 e 2ct ' s g{t/T) 2 dt 
s Jo 

increases (not necessarily strictly) as r increases, and is strictly positive for all r > unless y = 0, 
which can happen only if c = zq and r(t) = — c 2 /s. Integrating p.ip by parts, 

I c (r) = -y(0) - - / y(t)g(t/T)g(t/T) dt - - e~ 2ct ' s g(t/r) 2 dt - R c (r) . (3.2) 
T Jo s io 

As y is continuous, both e ct ^ s y(t)g(t/T) and e~ ct ^ s g(t/T) are square-integrable, and the Cauchy- 
Schwarz inequality gives the estimate 

1/2 /poo \ 1/2 

y{t)g(t/T)g(t/T)dt 



<(- s f y 2 e 2c ^g(t/r) 2 dt^J (s J™ e-**/'g(t/r)* dt^J 
< C^T~Rjr)e~ CT ' s , 



where C > is defined by 

C 2 =s I g(t) 2 dt 



and we have used the fact that git/r) = in [0,r]. We may also estimate 

^(l-e- 2c ^)<- j°° e-^g{t/rfdt<^. 
2 s J 2 

Putting these estimates in Q3.2p . we now have 

Ic(r) < -y(0) - |(1 - e- 2 -/ s ) - i? c (r) + 2C7e— /s ^5 



G 



and, estimating 



V T 

by the AM-GM inequality, we arrive at 



r R (t) p -2ct/s 

Io(t) < -1/(0) - £(1 - e" 2 ^) - ^ + 2 ^ 



2 V '2 t 

If either (a) r(i) is nonconstant, or (b) r(t) = r but zq ^ c or c 2 /s + r ^ 0, then ?/ cannot be 
identically zero. It follows that R c {t) is strictly positive and increasing on (0,oo). Thus there 
exists To > such that 

7 c (r) < -2/(0) - | = | - z 
for all r > To, which proves the required result. □ 

4 QEI-inspired hypotheses 

Now let us consider a QEI-inspired condition: we suppose that the forcing function r(t) is con- 
strained by inequalities of the form 

(r(t)~r (t))f(t) 2 dt>- HI/HI 2 (4.1) 

for all real- valued / G C XI (M), where ro(i) is a fixed continuous real- valued function and |||-||| is 
any Sobolev (semi-)norm of the form 

iii/iii 2 = X>ii/ (£) ii 2 , 

where Qi £ [0, oo) are constants of the appropriate dimensions, L £ N and || • || is the L -norm. 

In the analogy with QEIs, ro would be obtained from a reference state. As already mentioned, 
this constraint is inspired by, rather than exactly coinciding with, the known QEI results. See 
Section [5] for more discussion. 

Let g be as above and c > 0. Choose a smooth real- valued function h £ C°°(R), such that 
supp/i C [— To,oo) and h(t) — e~ ct l s on [0, oo). Then for each r > 0, 



Mt) = 



e~ ct/s 9(t/r) t > 
h(t) t < 



defines a test- function f T £ C X3 (M.). If c > then the dominated convergence theorem entails that 
as t — y oo and we have 



/0 />oo 
h(t) 2 (r(t) -r (t))dt- \\\h\\\ + liminf / e- 2ct/s r (t)g(t/T) 2 dt . 
-TO T^ + OO J Q 

The same holds if c — 0, provided that the coefficient Qo vanishes in the definition of |||-|||. Note 
that \\\h\\\ can be written 



The following result is now immediate. 



= E^fkc/ S ) 2 ^+ f ifcW(t)i s 

«_n \ Z •'-To 
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Theorem 4.1. Suppose r(t) obeys (|4.1[) and is nonconstant. Suppose there exist c > 0, tq > 
and h £ C°°(R) obeying supp/i C [— tq,oo) and /i(i) = e~ ct / s on [0, oo), /or which 



z - £+liminf / e- 2ct/s r (t)g(t/T) 2 dt> I h(tf Mt) - r (t)) dt + 

TTien (|2.4p /ias no solution on all of [0, oo). If r(t) = r obeys (|4.ip . £/ie same conclusion holds 
provided c ^ zq or c 2 / s + r ^ 0. // Qq = 0, the hypothesis on c can be relaxed to c > 0. 

Proof. By Lemma T3.ll if there is a solution to (12.41) on [0,oo) then 
> z a — — + liminf I c (t) 

2 T — >QO 

> z - - +liminf / e~ 2ct/s r (t)g(t/T) 2 dt- / h(t) 2 (r(t) - r (t)) dt - |||/i||| 2 

2 r^+oojp J_ Tq 

> 

by hypothesis, which is a contradiction. □ 
The above theorem tells us that a focal point will occur if there is sufficient initial contraction. 
Note that the dependence on the actual matter distribution r(t) is purely through its values in 
[— To,0], i.e., before the contraction zo is measured. At first sight it is strange that positive values 
of r{t) — rg require a larger value of zq: the reason is simply that the averaged nature of the 
energy condition means that large positive energies in the present allow large negative values in 
the future. The magnitude of initial contraction required is otherwise determined by the reference 
ro and the function h. 

In some QEIs, the support of the sampling function is constrained, for example, to be small 
relative to local curvature length scales. For simplicity, let us suppose that there is some ro such 
that (|4.1[) is valid for all / with support in an interval of length at most 2tq. (This would correspond 
to conditions in which curvatures remain bounded; if they do not, this might be regarded as an 
indication of singular behaviour, albeit not in the sense of geodesic incompleteness.) We also 
restrict to a particular seminorm which would be most relevant to QEIs in four dimensions. 

To discuss averages over longer timescales we will use a partition of unity. To this end, we 
will say that ip £ Cq°(R) is a bump function if it is nonnegative, supported in [0,3/2] and obeys 
V»(3/4 + x) = -0(3/4 - x), ip{x) 2 + 0(1/2 - x) 2 = 1 for x £ [0, 3/4]; in particular, this gives -0 = 1 
on [1/2, 1]. Then the functions ip(x ~ n) 2 (n £ Z) form a partition of unity on R. Note that this 
partition involves the squares of the underlying functions. The main work is to show that this 
gives usable bounds on the exponentially damped averages appearing in Lemma 13. II 



Proposition 4.2. Let ro be a fixed continuous function and suppose that r{t) is nonconstant and 
obeys 

/>oo 

(r(t)-r Q (t))f(t) 2 dt>-Q\\r\\ 2 (4.2) 



for any f £ Co°((0, oo)) supported in an interval of length at most 2tq > 0. Let ip be a bump 
function. Then for any c > we have 

/■oo poo rO 

liminf / e- 2ct/s r(t)g(t/T) 2 dt > liminf / e - 2ct/s r {t)g(t/T) 2 dt - (r(t) - r Q {t))h{t) 2 dt 

-«(i*iS(;)'+*<i^)- <«» 

for all real-valued h £ Cg°((-T ,Tb/2)) with h{t) = e" c */ s '0(l/2 - t/r ) on t > 0. The last term 
may be replaced by 

r 



2 (l-e- 2cr o/ s ) 



S 



We defer the proof to the end of this section. The following consequence is immediate, by the 
same reasoning as in Theorem 14. II 

Theorem 4.3. Let ro be a fixed continuous function and ip be a bump function. Suppose r(t) is 
nonconstant and obeys (|4.2[) for any f G Cq°((0,oo)) supported in an interval of length at most 
2tq > 0. If, for some c > and h obeying the conditions of Prop, we have 

z *>\- (RHS of O) (4.4) 

then (|2.4[) has no solution on all of [0,oo). If r(t) = r obeys (|4.2p . the same conclusion holds 
provided c ^ zq or c 2 / s + r ^ 0. 

Remark. It is possible to find bump functions ip with ||V'"|| 2 ~ 331. To see this, consider the ansatz 

{sin(0(x)) x < 1/4 

cos(0(l/2-x)) 1/4 < x < 3/4 
tp(3/2 - x) x> 3/4 

for 9 G C°°(M), 9 = on Mr, 9 > and 0(1/4) = tt/4, and 9^ 2k \l/A) = for all k G N, which 
ensure that rp is smooth at x = 1/4. After a calculation, one finds that 

U"\\ 2 = 2 I''" (9'{xf + 9"{xf) dt, 
Jo 

which corresponds to the Euler-Lagrange equation 

9""{x)-69'(x) 2 9"(x) =0 (4.5) 

together with boundary conditions 9{0) = 0'(O) = 0, 6»(l/4) = tt/4, 6»"(l/4) = 0. Numerical 
solution of this ODE in Maple 14 gives the value IIV 1 "!! 2 = 330.97. This does not extend to give 
a smooth solution to our original problem as 0"(O+) = 44.56; however, by considering a mollified 
function 6 e (x) = 6(x)f(x/e) with / G C°°(M), f(x) = on (— oo,0], f(x) = 1 on [l,oo), one may 
show that ||Ve II ~> IK II as e -> 0+ [one uses estimates 0(x) = 0(x 2 ), 0'(a;) = O(x), 0"(x) =0(1) 
near x = 0]. 

Proof of Proposition \4-S\ The idea of the proof is to use a partition of unity to decompose the 
integral 

POO 

Jc(t)= / e~ 2ct / s (r(t)-r (t))g(t/T) 2 dt, 
Jo 

thereby obtaining a sum of integrals, each of which can be bounded below using (|4.2|) (only 
finitely many terms appear at any fixed value of r). The sum of the resulting lower bounds must 
be controlled as r oo. 

We start from the identity 

For real- valued ip G Cq°((0, oo)), we may integrate by parts and discard a nonpositive term to find 
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If, additionally, the support of ip is contained in an interval of length at most 2 To, we may 
apply (|4.2j) with f(t) — e~ ct / s g{t/T)ip(t) and g chosen as above, to give 



oo 

e 



- 2 ^ s (r(t)-r Q (t))g(t/r) 2 ^t) 2 dt 

4 / ,;J 



> -Q^°V 2ct / s i^git/rf^tf + (J^gWrMt)) j dt, (4.8) 

where we have estimated ||/"|| 2 using (|4.7[) , with y> replaced by g(t/T)tp(t). This estimate will be 
used for each term appearing in the partition of unity, to which we now turn. 

Let tp(t) = ^>(£/to), where ip is the chosen bump function. Defining tp n (t) = tp(t — nto), the 
(fin form a partition of unity so that J2r7=o fn{t) 2 = 1 on t > tq/2 and at most two tp n {t) are 
nonzero at each t; each ip n has support diameter of 3tq/2. We may now decompose the integral 
J c (t) using this partition, estimating each term using (|4.8[) . This gives 

/■To/2 °° />oo 

Mr) = / e- 2ct ' s {r{t) - r (t))g(t/r) 2 (1 - ^) 2 ) dt + £ / e- 2ct ' s (r{t) - r (t))g(t/r) 2 <p n (t) 2 dt 
J n = J0 
rro/2 oo 

e-**/'(r(t) - r (t))g(t/T) 2 (l - V {tf) dt-Q^ S n (r) , 

n=0 



> 



where 



Our task is now to control the sum X)^°=o ^«( r ) m ^ ne Ihrfit r — > oo. Now, for r > (n + 3/2)r , 
g(t/r) = 1 on the support of <p n , and we obtain 

4 \ 

tit 



5n(r) = jf e- 2 ^ (t) 2 + ^(t) 2 ) 

< (£) 4 H e - 2ct / s (p n (t) 2 dt + e- 2cnT «/ s \\<p"\\ 2 , 



where the exponential factor in the second term arises on account of the support properties of tp n 
and we use ||<p^|| = Hv 5 "!!- Provided the limit and sum can be exchanged (which will be justified 
below) we then have 



e -2cnT /s 



r lim £S n (r)<£Q) / e- 2 «l^n{t) 2 dt + \W\\ 2 Y, 

The first sum can be evaluated using the partition of unity property and the second is elementary, 
so 

4 r T °/ 2 , /r\i r°° , \\w"\\ 2 



Urn J2Sn(r)<(-) / e ~ 2 ^Mt) 2 -l)dt+^) / 



n=0 

„"l|2 



1 _ g— 2cto/s 



2 Vs/ 1 - e - 2cT <i/ s ' 

where we discard the nonpositive contribution from <y9 2 — 1 in the last step. 
Assembling the results so far, we have 

r°l 2 i /I /r\3 lh//'ll 2 

liminfj c (r)> / e~ 2ct/s {r{t) - r (t))(l - ip(t) 2 )dt - Q ( ~ (-) + W 



2ct /s) 
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Taking h as specified in the hypotheses, we note that h(t) 2 = e 2ct / s {\ — ip(t) 2 ) on t > 0. Applying 
Eq. (|4.2[) with /i in place of /, we obtain 



liminf Ur) > - j_Jr(t) - Mt))h(t) 2 dt-Q (||^'|| 2 + \ + 

/" OO 

liminf J c (r) > liminf J c (r) + liminf / e" 2c * /s r (i) eft , 

r— 5-+oo t— >+oo r— >+oo Jq 



Finally, because 



we obtain the the required result (|4.3p . The alternative form of the bound arises if we estimate 
\\h"\\ 2 by using (14.61) applied to ip(t) = tp(l/2 — t/r ) and integrating by parts on [0, 00) to obtain 
(|4.7[) (the boundary terms at t = cancel). Using the properties of ip we find 



\h"\\ 2 < 



\ h " { t)\ 2 d t + \^\'w> 



which yields the alternative bound. 

It remains to justify the exchange of limit and sum used above, for which we use a dominated 
convergence argument. We split S n (r) into the sum of two terms. For the first, we note that 

00 

e - 2ct/s g(t/T) 2 (p n (t) 2 dt < e- 2cnT °/ s M 2 

because aupp ip n C [nrg, 00). The second is regarded as the square of an L 2 norm, which can be 
estimated using the Leibniz rule and triangle inequalities as 



" ' ' ^f(*/r)Vn(t)) dt 



1/2 



where C n , C' n and C[[ are the L 2 norms of e~ ct / s g"{t/T)ip n (t), 2 e - ct / s g / (t/T)(p' n (t) and e~ ct / s g{t/T)if 
respectively. Setting C — max{l, 4||g'||^ c , H^'H^} and again using the support properties of ip n , 
we have 



r e ~ 2ct/s (^ 5(t/r) ^ ( ^) * - c ^~ 2cnT ° ,s 



l^'ll , IIV'II 



TTQ 



Accordingly, there exists K > such that < 5„(r) < Ke 2cnT <>/ s for all n e N and r > r , 
where K is independent of r and n. As this upper bound is evidently summable, the dominated 
convergence theorem permits us to exchange the limit and sum as claimed. □ 



5 Singularity theorems 

The results from the previous sections can be used to obtain refinements of some of the classical 
singularity theorems of Hawking and Penrose [H] that allow global violations of the classical 
energy conditions. 

We first consider the cosmological setting. 

Theorem 5.1. Let M be a globally hyperbolic spacetime of dimension n > 2, and let S be a 
smooth compact spacelike Cauchy surface for M . Suppose along each future complete unit speed 
timelike geodesic 7 : [0, 00) —¥ M issuing orthogonally from S, there exists c > such that, 

liminf / e - 2ctl{n - l) r{t) dt > 9(p) + - , (5.1) 

where r(t) := Ric (7'(t),V(£)) = Rab"f' a J lb (t) and 9(p) is the expansion (i.e., mean curvature) of 
S at p = 7(0). Then M is future timelike geodesically incomplete. 
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Remarks. 



1. As usual, if one assumes that the Einstein equations, 

Rab - ^Rg a b = 8irT ab (5.2) 

hold then the Ricci curvature term can be directly related to the energy-momentum tensor 
T a b, and, hence, one should view (|5.1[) as an energy condition on spacetimc. 

2. Note that the energy condition (|5.ip is satisfied provided along each such timelike geodesic 
7, the condition, 



liminf / r(t) dt > 9(p) 



(corresponding to c = 0) holds. This condition emphasizes the fact that if S is mean 
contracting, i.e., if 9 is negative on S, then (|5.1I) can hold even if r(t) is everywhere negative. 

3. Specializing further, if the strong energy condition holds, i.e., if Ric (X,X) = R ab X a X b > 
for all timelike vectors X, and if S is everywhere mean contracting then (15.11) holds, and 
we essentially recover Hawking's cosmological singularity theorem; cf., Theorem 4 in [211 p 
272]. (For convenience we have stated Theorem 15. II as a future singularity result rather than 
a past singularity result.) 

Proof of Theorem \5. 11 We construct an S'-ray 7, i.e., a future inextendible timelike geodesic 7 
emanating from a point on S that realizes the Lorentzian distance to S from each of its points, 
as follows. Choose a sequence of points q n which extends arbitrarily far into the future of So 
By properties of Cauchy surfaces, there exists a timelike geodesic segment 7„ from p n G S to q n 
that realizes the Lorentzian distance from S to q n . Since S is compact, by taking a subsequence 
if necessary, we may assume that the sequence p n converges to a point p G S. Let 7 : [0, a) — > M, 
a G (0,oo], be the future inextendible unit speed timelike geodesic issuing orthogonally from 
p G S. The maximality of the 7„'s guarantees that 7 is an S-ray (see, for example, the proof of 
the Sublemma in [25]). 

Let p : J + {S) — >• M be the Lorentzian distance function from S, 

p(x) = d(S,x) = supd(y,x) . 

yes 

By global hyperbolicity, p is continuous on J + (S), and will be smooth up to the focal cut locus 
of S (see [26 ). Since 7 is an S'-ray, there are no focal points, or focal cut points, to S along 7, 
which, by the lower semi-continuity of the S'-distance-to-cut locus function [2 6) , ensures that p is 
smooth on a neighbourhood U of 7. On this neighbourhood u := — Vp is a smooth future directed, 
geodesic, hypersurface orthogonal unit timelike vector field such that u = 7' along 7. Consider 
the expansion scalar 9 — divu. Along 7, 6 = 9(t), t G [0, a), obeys Raychaudhuri's equation (for 
an irrotational timelike congruence) |21j . 

§ = -Ric( 7 ',7')-2a 2 -^-fl 2 , (5.3) 
dt n — 1 

where a is the shear scalar. 

We now observe that 7 is necessarily future incomplete. For suppose that 7 is future complete 
(i.e., suppose a = 00). Setting z — —9, r = Ric (7', 7') + 2er 2 , s — n — 1 and z(0) = —9(p) in 
(I2.4|) . we see that Lemma [2.31 together with the energy condition (|5.1j) . implies that (|5.3[) . with 
9(0) = 9(p), has no solution on [0, 00), which is a contradiction. □ 

Remarks. 



4 More precisely, if h is a complete background Riemannian metric on M, choose q n in J + (S) so that the 
fc-distance from S to q„ tends to infinity as n — > 00. 
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1. Although slightly more complicated to state, Lemma I3TT1 and Theorems 14. II and 14.31 provide 
alternative energy conditions that yield singularity theorems similar to Theorem 15. II 

2. The assumption of global hyperbolicity in Theorem 15.11 can be relaxed. It is sufficient to 
assume that M admits a smooth compact acausal spacelike hypersurface S. In this case 
one can construct an 5-ray 7 contained in the future domain of dependence D + (S) (see 
especially [25, Main Lemma]), and use the fact that S is a Cauchy surface for the total 
domain of dependence D(S). 

3. We see from the proof that it would be sufficient for (|5.ip to hold on S'-rays. 

We now consider an extension of the Penrose singularity Theorem [2TJ Theorem 1] to the 
energy conditions being considered here. Let S be a codimension two compact acausal spacelike 
submanifold in a spacetime M of dimension n > 3. Under suitable orientation assumptions, 
there exist two smooth nonvanishing independent null normal vector fields, £ + and £-, along £, 
corresponding to outgoing and ingoing light-rays emanating from E. The null expansion scalars 
9+ and 9- on E are obtained by taking the divergence of £+ and £-, respectively, along £, and 
measure the instantaneous divergence of the outgoing and ingoing light rays emanating from E. 
In a strong gravitational field both 9+ and 9- can be negative, in which case E is called a trapped 
surface. According to the Penrose singularity theorem, if E is a trapped surface in a spacetime M 
having a noncompact Cauchy surface and satisfying the null energy condition, then M is future 
null geodesically incomplete. 

The following generalizes the Penrose singularity theorem (see also [UJ Theorem 2] ) . 

Theorem 5.2. Let M be a spacetime of dimension n > 3 with a noncompact Cauchy surface 
S. Let E be a smooth compact acausal spacelike submanifold of M of codimension two, with null 
expansion scalars 9± associated to the future directed null normal vector fields £± . Suppose along 
each future complete affinely parameterized null geodesic r\ : [0, 00) — > M , issuing orthogonally 
from E with initial tangent £± , there exists c > such that, 



where p — f](0) and r(t) := Ric (r]'(t), r}'{t)) — R a br} la rj' h {t). Then M is future null geodesically 
incomplete. 

Remarks: If E is a trapped surface and M obeys the null energy condition, then (|5.4[) is satisfied 
for sufficiently small c and we recover the Penrose singularity theorem. The case in which (|5.4p 
holds with c = corresponds to the codimension two case of Theorem 2 in [27]. As in Theorem 
15.11 Lemma 13.11 and Theorems 14.11 and 14.31 may be used to provide alternative energy conditions 
that yield singularity theorems similar to Theorem l5.2l Finally note that if the Einstein equations 
(|5.2[) hold, then the Ricci curvature term Ric (rj'(t), n'(t)) can be replaced by 8irT ab r]' a r]' b {t). 

Proof of Theorem \5.SX We shall be brief, as the general structure of the proof is similar to the 
proof of the Penrose singularity theorem. See |21j for further details and relevant results from 
causal theory. 

The achronal boundary 9J + (E) is a C° achronal hypersurface ruled by null geodesies, i.e., for 
each point q € <9J + (E) \ E there exists a future directed null geodesic segment from a point in E 
to q which is entirely contained in <9J + (E). Such null geodesies are called the null generators of 



If 9J + (E) were compact, then flowing along the integral curves of a timelike vector field on 
M would establish a homeomorphism between 9J + (E) and the Cauchy surface S, which would 
contradict the noncompactness of S. Thus we may assume 9J + (E) is noncompact. Since <9J + (E) 
is closed but noncompact we can find a sequence of points q n in dJ + (E) whose /i-distance to E 
(where h is a complete background Riemannnian metric on M) tends to infinity as n — > 00. Let r\ n 
be a null geodesic generator of 9J + (E) fromp„ e E to q n ; i] n must meet E orthogonally, otherwise 




(5.4) 



&/+(£). 
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the achronality of <9J + (E) would be violated. By passing to a subsequence if necessary, we may 
assume without loss of generality that p n — > p G E and that the initial tangents tf n (0) —> i+(p)- 

Let r\ : [0,a) —¥ M, a £ (0, oo], be the affinely parameterized future inextendible null geodesic 
emanating from p with initial tangent £+(p); rj is contained in <9J + (E) since each rj n is. Since 
SJ + (E) is achronal there can be no null focal points, or null focal cut points, to E along rj. It 
follows [26] that rj is contained in a smooth null hypersurface H C 9J + (£) generated by null 
geodesies emanating from £ near p with initial tangents given by £+. Let K be a smooth null 
vector field on H such that K — rj' along 77, and let be the null expansion of H with respect to K. 
Along 77, 9 — 9(t), t g [0,a), obeys Raychaudhuri's equation (for an irrotational null congruence) 

ED, 

§ = -mc{r 1 '^)-2^-^-e\ (5.5) 
at n — I 



Now one can argue just as in the proof of Theorem 15. II If 77 were future complete, so that a — 00, 
then Lemma \2. 31 and (|5.1I) would imply that (|5.5p . with 6>(0) = 0+(p), has no solution on [0,oo). 
Hence, rj must be future incomplete. □ 



6 Applications and discussion 

We begin with two applications to the Einstein equations coupled to a real scalar field. The 
minimally coupled Einstein-Klein-Gordon system in n > 2 spacetime dimensions is given by 

Rab - \Rgab = 8^C in , = Va4>^b<t> ~ ^at(V^V^ + m 2 <}?) . 

If 7 is a unit speed timelike geodesic then 

r(i) = Ric( 7 ',7') = 8^ ((V y 0) 2 - -^-yf) > 

from which we can see easily that (see, e.g., [28] and [21] p. 95]) the Klein-Gordon field fails to 
obey the strong energy condition and hence the standard hypotheses of the Hawking singularity 
theorem. However, if 4> remains bounded along 7, with \<j>\ < <fi*, say, then we have 

- - + [ T e-^/f^Vft) dt>--- — 
2 Jo y > ~ 2 2c 

for all c,T > 0, where K = m<f)* ^/87r(n — l)/(n — 2). The right-hand side has a maximum for 
c = K , and we may conclude from Theorem 15 . 1 1 that if 8 < —K everywhere on 5* then M is future 
timelike geodesically incomplete. 

In fact, the same result follows from Theorem 4 of |15] (in the n = 4 case). However, we may 
also obtain more general results, which both allow for exponential growth of <f> in the proper time 
along these geodesies and provide finer detail in the case of exponential decay. For example, if S 
is a smooth compact spacelikc Cauchy surface and the scalar field obeys a bound 

1^0)1 < (^e^/t"- 1 ) 

for p € J + (S), where p(p) is the Lorentzian distance from p to S and a £ R is constant, then we 
obtain a bound 

-- + [ T e - 2ct /(«-i) r (t) dt > -- - R2 , 
2 i W " 2 2(c - a) 

for any c > a and all T > along S'-rays, for which p("f(t)) — t. Optimising over c > as before, 
and using Theorem 15. II (and the third remark following its proof) we then have that M is future 
timelike geodesically incomplete if 




-a/2-K a>-K 
K 2 /{2a) a< -K 
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on S. The case a = corresponds to the result above, while the other cases indicate the additional 
power of our technique. 

A drawback of these results is the need to invoke the supremum (/>* and (if needed) the con- 
stant a. In principle they are implicit in the Cauchy data for the Einstein-Klein-Gordon sys- 
tem on S. One can also read this result in the following way (with a = for simplicity): if 
S is mean contracting then either M is future timelike geodesically incomplete or <f> exceeds 
TO_1 \/ (87r) _1 (n — 2)/(n — 1) infs \0\ in magnitude somewhere to the future of S. 

We now turn to the nonminimally coupled field of mass m > and coupling £, with stress- 
energy tensor 

Tab = T$» + £ (g ab D - V a V b + {R ab - ^Rg ab )\ <\> 2 ■ 

The presence of the Einstein tensor in this expression complicates the initial value problem for the 
Einstein equations with this source (see [29] for existence and uniqueness results). In addition, 
the second derivatives allow violations of the NEC as well as the SEC, which can be exploited 
to construct nonsingular cosmological solutions [30. However, for £ g [0,1/4] (including the 
conformal coupling £ = \{n— 2)/(n — 1)) the theory obeys the inequality 

J WVV(A) 2 d\ > -2£ J {/' (A) 2 - \mc (V, 7 ')/(A) 2 - ( \ dj Ri 2 ) dX 

for any affinely parameterized causal geodesic 7 and smooth compactly supported, real- valued / 
Q311 Theorem II. 1], modulo change in conventions). If we assume that the field magnitude remains 
bounded, \</>\ < </>* this entails (for the case of null 7) that 



Ric (7', 7 ')(1 - 8^ 2 )/(A) 2 dX > -16^ 2 ||/'ll 

for all / € Cq°(R). Provided that 0» is strictly less than the critical value (87r^) -1 / 2 we may 
absorb a factor into /, thus obtaining 

Ric(7',7')/(A) 2 dA>-167r^ 2 [°° (±-J2L=\ dX 



>-o(ii/'ii 2 +Q 2 n/ir) 



for all real- valued / € Cq°(R), where 



Q _ 32^0 2 - = 8tt^*^ 



1 - Stt^^ 1 - 8tt^ 2 

and is an upper bound on \(j)'\. (We have used the simple estimate ||(/g)'|| 2 < 2(||/|| 2 ||.g'|| 2 K; + 
||/' IPNIL) for any / € C* °°(R), g € C°°(R).) This is a bound of the form gT]) with r = and 
r = Ric (7', 7')- To proceed, we need the following lemma. 

Lemma 6.1. Let t > 0, s > and choose K > so that 

K 2 >Q 2 + Q- 1 Ric (7', 7') 
on (— ro,0]. For any e > 0, £/iere exists c > and ft, € C°°(R) supp/i C [—To, 00), h(t) = 



/■0 1 

/ h 2 Ric (7', 7') dA < Qy / Qs + Q 2 /2 + -QK coth Kt 

J -Tn ^ 
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Proof. We take c = Qsy/2Q / (Q + 2s). Then the left-hand side becomes 



Qy/Qs + Q 2 /2 + / (Oh 12 + [QQ 2 + Ric (7', i)]h 2 \ dX 
so it is sufficient to show that 

r° k 

inf / (ti 2 + K 2 h 2 ) d\<— cothifro , (6.1) 

with the infimum taken over the class of h specified in the hypotheses. Treating this as a varia- 
tional problem, the Euler-Lagrange equation is h" = K 2 h and applying the boundary conditions 
h{~T ) — 0, h(0) — 1, the solution is h(X) = (sinhifro)" 1 sinhi4T(A + r ), giving equality in (|6.1[) . 
This can be approximated arbitrarily well within the given class of h. □ 
Accordingly, using Theorem 14.11 in place of Lemma 12.31 we obtain an analogue of Theorem 15.21 
if the hypothesis (|5.4[) is replaced by the requirement 

0±{p) < ~Q^(n-2)Q + Q 2 /2- X -QK coth Kt q , 

with K computed as above for the extension of the null geodesic r\ to (— tq,0], for some To > 0. 

Finally, we briefly discuss the potential for applications involving quantum fields. The QEIs 
established for many free fields are weakened versions of the weak and dominant energy conditions, 
while our analogue of the Hawking singularity theorem is based on a weakened strong energy 
condition. In order to prove a Hawking result from QEI hypotheses, it would be necessary to add 
conditions on the trace of the stress-energy tensor, much as in the case of the classical scalar field 
discussed above. In addition, it is necessary to find estimates on the timescales for averaging over 
which the curved spacetime QEIs are well-approximated by Minkowski space results. Nonetheless, 
there seems a good prospect of obtaining results along these lines, and our results provide a proof 
of principle for the idea that energy conditions based on local averages can be used to deduce 
singularity theorems. Turning to the Penrose-type results, a significant problem at present is that 
no locally averaged energy inequality along individual null geodesies is known, and direct analogues 
of the results for timelike averaging cannot hold in dimensions higher than 2 [35]. Here, it seems 
that the best approach is to consider averages with a degree of transverse smearing, whereupon 
QEIs can be proved [32]. It is worth recalling that even the ANEC condition is problematic for 
averaging along a complete null geodesic: the real scalar field violates ANEC in general spacetimes 
(see [33J and references therein) although it holds in Minkowski space [T3J and in general spacetimes 
along complete achronal geodesies with a tubular Minkowskian neighbourhood [31]. Moreover, this 
problem persists even with some transverse averaging schemes }35). 

Finally, it should be noted that a full treatment of singularity theorems in the context of 
quantized matter would require (at least) a semiclassical analysis that takes backreaction into 
account in a dynamical fashion, bringing significant technical challenges - see |36[ 137) for recent 
results in simple cosmological models and further references, and |38j for positive results on ANEC 
with transverse smearing in this context 
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